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The rapid advance in experiments with rare isotopes has opened up several new frontiers 
in nuclear science [0-|^. It is now possible to study in detail structures of nuclei in many 
unexplored regions of the periodic chart and novel properties of highly isospin asymmetric 
nuclear matter. Prospects for new physics in nuclear reactions induced by especially exotic 
neutron-rich nuclei have generated much interest in the nuclear science community. In par- 
ticular, fast fragmentation beams at the planned Rare Isotope Accelerator (RIA) will provide 
a unique opportunity to explore features of nuclear matter at extreme isospin asymmetries 
towards the pure neutron matter 0]. The isospin-dependence of the nuclear equation of 
state (EOS) is among the most important but very poorly known properties of neutron-rich 
matter Pj. It is relevant to Type II supernova explosions, to neutron-star mergers, and to 
the stability of neutron stars. It also determines the proton fraction and electron chemical 
potential of neutron stars at (3 equilibrium. These quantities consequently determines the 
cooling rate and neutrino emission flux of protoneutron stars and the possibility of kaon 
condensation in dense stellar matter [p|-p!^. At present, nuclear theories predict vastly dif- 



ferent isospin-dependence of the nuclear EOS depending on both the many-body techniques 
and the bare two-body and/or three-body interactions employed. In this Letter, we show 
that the neutron-proton differential flow is a very useful probe of the isospin-dependence of 
the nuclear EOS. This approach utilizes constructively both the isospin fractionation and 
the nuclear collective as well as their sensitivities to the isospin-dependence of the nuclear 
EOS. It also avoids effectively uncertainties associated with other dynamical ingredients of 
heavy-ion reactions at intermediate energies. The experimental measurement of neutron- 
proton differential flow provides thus a novel means for determining the isospin-dependence 
of the nuclear EOS. 



Various theoretical studies (e.g., [p!3| , p!4| ) have shown that the energy per nucleon e(p, 5) 
in nuclear matter of density p and isospin asymmetry parameter 5 = (p„ — Pp)/{pn + Pp) can 
be approximated very well by a parabolic function 

e(p,<5) = e(p,0) + 5(p)5^ (1) 



where e(p, 0) is the EOS of symmetric nuclear matter. In the above S{p) is the symmetry 
energy 

5(p) = (2^/^-l)^i?V/^ + 5pot(p), (2) 

where Ep is the Fermi energy at normal nuclear matter density po and 5'pot(p) is the poten- 
tial contribution to the symmetry energy. Vastly different density-dependences of Spot{p) 
given by various many-body theories have led to rather divergent predictions on the isospin- 
dependence of the nuclear EOS. Here we adopt the parameterization of S'pot(p) from the 
study of neutron stars (e.g., [p|, |l5| , p^ ) 



SMp) = {Sipo) - (2^/3 - l)^i^°) Fiu), (3) 

where u = p/po is the reduced density and S'(po) is the symmetry energy at po. The 
latter is known to be in the range of about 27-36 MeV |[T8|-pO|. We consider two typi- 



cal parameterizations of Spotip) with F{u) = 2m^/(1 + u) and F{u) = v}!'^ , respectively. 
The isospin-dependence of the nuclear EOS can be characterized by the curvature of the 
symmetry energy at po 



2 5'5(P) 
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For the two parameterizations of 5'pot(p) considered here the K^ym parameter is 61 MeV and 
—69 MeV, respectively. These values are well within the wide range of K^y^ from about 
—400 MeV to -1-466 MeV predicted by many-body theories (e.g. HJT3[]). Available data of 
giant monopole resonances does not give a stringent constraint on the -ft"sym parameter either 
p^ . The corresponding single particle symmetry potential Vl^-^ is 

2 2 

and 

K:/^ = ^le^u'lH - \eyH\ (6) 



for Ksym = +61 MeV and —69 MeV, respectively. In the above Ca is a constant of 

ea = S{po) - (2^/3 _ 1) 1^0 ^ ig MeV. (7) 

The + and (— ) sign is for neutrons and protons, respectively. For small isospin asymmetries 
and densities near po the above symmetry potentials reduce to the well-known Lane potential 
which varies linearly with 6 [p4 |. 



Our study is based on the isospin-dependent Boltzmann-Uehling-Uhlenbeck (IBUU) 
transport model (e.g., |]2llj2^ . In this model protons and neutrons are initialized according 



to their density distributions predicted by the relativistic mean-field (RMF) theory p3[ . The 



isospin-dependent reaction dynamics is included through isospin-dependent nucleon-nucleon 
scatterings and Pauli blockings, the symmetry potential V^(^^ and the Coulomb potential 
VJ' for protons. A Skyrme-type parameterization is used for the isoscalar potential. For a 
recent review of the IBUU model, we refer the reader to ref. 0]. 

To establish the essence and to reveal the advantage of analyzing the neutron-proton dif- 
ferential flow, we first study individually the isospin fractionation and the isospin-dependence 
of nucleon collective flow. These phenomena in nuclear reactions induced by neutron-rich 
nuclei are very interesting in their own rights. Here we concentrate on studying their sen- 
sitivities to the isospin-dependence of the nuclear EOS only. The isospin fractionation is 
an unequal partitioning of the neutron to proton ratio N/Z of unstable asymmetric nuclear 
matter between low and high density regions. It is energetically favorable for the unstable 
asymmetric nuclear matter to separate into a neutron-rich low density phase and a neutron- 



poor high density one p5|-p9|]. This phenomenon has recently been confirmed in intermediate 
energy heavy- ion experiments [^0|,^. We quantify the degree of isospin fractionation by 
using the ratio of (N/Z) free to {N/Z)hound- Here (N/Z) free and {N/Z)bound are the saturated 
neutron to proton ratios of nucleons with local densities less (free) and higher (bound) than 
Pc = l/Spo, respectively. The density cut-off pc has often been used in transport models for 
identifying approximately free nucleons quickly. The specific value of pc used here does not 
affect our conclusions in this work. Shown in Fig. 1 is the degree of isospin fractionation 



as a function of the neutron to proton ratio {N/Z)sys of the reaction system (left), impact 
parameter (middle) and beam energy (right), respectively, for reactions between several Sn 
isotopes. It is seen that the degree of isospin fractionation increases with both the {N/Z)sys 
and the impact parameter, but decreases with the beam energy. It is very interesting to 
see that the degree of isospin fractionation is rather sensitive to the Kgym parameter. This 
sensitivity is very strong noticing that the change of Ksym made here is only about 1/7 of its 
variation predicted by the many-body theories. The origin of isospin fractionation and its 
dependence on the Kgym parameter can be easily understood from the density dependence 
of the symmetry energy. Since the repulsive symmetry potential for neutrons increases with 
density more neutrons will be repelled from high density regions to low density regions. 
While the opposite is true for protons because of their attractive symmetry potentials. Fur- 
thermore, the magnitude of symmetry potentials is higher for Ksym = —69 MeV than for 
Ksym = +61 MeV for densities less than about po- One thus expects to see a higher degree 
of isospin fractionation with Ksym = —69 MeV as shown here. 

The analysis of collective flow has been very useful in studying various properties of 
nuclear matter (e.g., P2|-p5|). How useful is it for studying the isospin-dependence of the 



nuclear EOS? To answer this question we study the average in-plane transverse momentum 
< Px/N > of free nucleons as a function of their reduced rapidity y/ycms- Shown in Fig. 2 
is such an analysis for the reaction of ^'^'^Sn +^^^ Sn at a beam energy of 50 MeV/nucleon 
and an impact parameter of 5 fm. To make the analysis more sensitive to the symmetry 
potential, the incident energy is chosen to be near the balance energy of the reaction system 
according to the systematics of balance energies [^]. Our calculations are performed with 



the Ksym parameter of +61 MeV (upper window) and —69 MeV (lower window), respectively. 
The difference between neutron- and proton-flow is found to depend strongly on the Ksym 
parameter. Thus, interesting information on the isospin-dependence of the nuclear EOS can 
also be obtained from analysing the collective flow of nucleons. In the case of Ksym = +61 
MeV the repulsive Coulomb potential dominates over the negative symmetry potential for 
generating the proton-flow. The net sum of Coulomb and symmetry potentials for protons 



is actually higher than the positive symmetry potential for neutrons. The proton-flow is 
thus higher than the neutron-flow. In particular, near the mid-rapidity, neutrons are still 
flowing to the negative direction because of the attractive mean filed acting on them while 
protons have already started flowing in the positive sense. By changing the Ksym parameter 
from +61 MeV to —69 MeV the magnitude of symmetry potentials is increased but the 
Coulomb potential remains the same. Thus the difference between neutron- and proton-flow 
gets much reduced as one expects and shown in the lower window of Fig. 2. 

It is well known that collective flow is also sensitive to both the curvature Kq of the 
isoscalar part of the nuclear EOS and the in-medium nucleon-nucleon cross sections (e.g., 
|39|-^). Both quantities are still associated with some uncertainties despite the remarkable 



progress made in the last two decades. Also the symmetry potential is normally rather small 

compared to the isoscalar one. Therefore, to study the isospin-dependence of the nuclear 

EOS one has to select delicate observables that minimizes influences of the isoscalar potential 

but maximizes effects of the symmetry potential. It would be ideal if these observables can 

also avoid effects of other dynamical ingredients of the reaction dynamics. We introduce 

here such an observable, i.e., the neutron-proton differential flow 

I N(y) 
Fn-p{y) = j;^ E P=^.^" (8) 

where N{y) is the total number of free nucleons at rapidity y and Tj = 1 and —1 for neutrons 
and protons, respectively. The neutron to proton ratio in N{y) is determined by the degree 
of isospin fractionation. In essence, this isospin-averaged nucleon flow combines construc- 
tively contributions of the symmetry potential to the collective flow of both neutrons and 
protons. Simultaneously, it reduces significantly influences of the isoscalar potential and the 
in-medium nucleon-nucleon cross sections as they both act identically on neutrons and pro- 
tons. Therefore, it uses efficiently the isospin effects shown in both the isospin fractionation 
and the nucleon collective flow. Shown in Fig. 3 is the neutron-proton differential flow as a 
function of the reduced rapidity for collisions between three Sn isotopes at a beam energy 
of 50 MeV/nucleon. A clear signal of the isospin-dependence of the nuclear EOS is seen at 

6 



both the target and projectile rapidities. The signal becomes stronger with the increasing 
isospin asymmetry of the reaction system as one expects. Moreover, we found that a 40% 
variation of the Kq parameter or a change by a factor of 2 of the in-medium nucleon-nucleon 
cross sections results in a less than 7% change in the values of Fn-p{y). The neutron-proton 
differential flow is thus a very useful probe of the isospin dependence of the nuclear EOS. 
We have also performed a study on the excitation function of the neutron-proton differen- 
tial flow from Ebeam/A = 30 MeV to 400 MeV for the Sn+Sn reactions at several impact 
parameters. The signal exists clearly in the whole energy range and is the strongest close 
to the lower energy end. For a given reaction system involving neutron-rich nuclei, what is 
the best beam energy to learn the most about the isospin-dependence of the nuclear EOS? 
We suggest that one selects the incident energy to be close to the balance energy of a re- 
action system with similar mass but along the (3 stability valley. It is well known that the 
transverse flow disappears at the balance energy where the attractive interactions balances 
the repulsive scatterings P5|-p^. At this energy the non-isospin related background is the 
smallest for studying the isospin-dependence of the nuclear EOS with the neutron-proton 
differential flow. It is useful to point out here that the well-established systematics of bal- 



ance energies |5^ can serve as a guide in selecting the best beam energies for the purposes 
discussed above. 

In conclusion, a novel means for determining the isospin-dependence of the nuclear equa- 
tion of state, i.e., the analysis of neutron-proton differential flow, is introduced for the flrst 
time. This approach uses constructively both the isospin fractionation and the nuclear col- 
lective flow as well as their sensitivities to the isospin-dependence of the nuclear EOS. It 
also avoids effectively uncertainties associated with other dynamical ingredients of heavy-ion 
reactions at intermediate energies. 
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FIG. 1. Isospin fractionation as a function of (N/Z) of the rection system (left), impact param- 
eter (middle) and beam energy (right) in Sn+Sn reactions 
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FIG. 2. The average transverse momentum per nucleon in the reaction plane for neutrons and 
protons as a function of reduced rapidity with the Kgym parameter of +61 MeV (upper window) 
and —69 MeV (lower window), respectively. 
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FIG. 3. Neutron-proton differential flow for the reaction of ^^'^Sn -j-^^^ g^ ^124 ^^ _j_i24 g^ ^^^ 
^'^'^Sn +132 g^ ^^ ^ beam energy of 50 MeV/nucleon, an impact parameter of 5 fm and the Kgy^ 
parameter of +61 MeV (filled circles) and —69 MeV (open circles), respectively. 
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